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whereas anotherLSBBEZI finds 9 in the range 0.26 - 0.31. 
Ground state methods involve computing the ground 
state with two different boundary conditions, and one 
of the main reasons for the large range of estimates for 9 
is that different groups made different assumptions about 
the optimal form of boundary condition changes, see, 
e.g., Ref. |lj. This difficulty is avoided in the finite- 
T Monte Carlo approach adopted here. We find that 
9 = 0.27 ± 0.01, which agrees with the results of the sec- 
ond group of ground state calculations mentioned above. 

In Sec. H we introduce the model and observables 
as well as details regarding the Monte Carlo technique 
used. Some results from finite-size scaling are discussed 
in Sec. [II. Our results in two and three dimensions are 
presented in Sees. |fV] and |v|, respectively, and the con- 



We present results from Monte Carlo simulations of the two- and three-dimensional gauge glass at 
low temperature using the parallel tempering Monte Carlo method. Our results in two dimensions 
strongly support the transition being at T c — 0. A finite-size scaling analysis, which works well only 
for the larger sizes and lower temperatures, gives the stiffness exponent 9 — —0.39 ± 0.03. In three 
dimensions we find 9 — 0.27±0.01, compatible with recent results from domain wall renormalization 
group studies. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

The gauge glass is a model which is often used to de- 
scribe the vortex glass transition in high-temperature su- 
perconductors. Consequently there has been a substan- 
tial amount of theoretical work attempting to understand 
its properties. Nonetheless, there are still areas of dis- 
agreement. For example, in two dimensions there is an 
ongoing controversy as to whether a spin-glass transition 
occurs at finite temperature or not. In addition, in three 
dimensions, although it is known that the system exhibits 
a finite-temperature transition, see, e.g., Ref. [|, there is 
no consensus on the value of the stiffness exponent. 

Recent results by means .of resistively-shunted junc- 
tion (RSJ) dynamics by Kimo claim evidence for a finite- 
temperature transition in two dimensions with a tran- 
sition temperature T c = 0.22. Choi and ParkEl find the 
same value for T c by studying the scaling of the spin-glass 
susceptibility via Monte Carlo simulations. These claims 
for a finite Tq-are in contrast to results by GranatoQ and 
Hyman et al.a who also use RSJ dynamics and find ev- 
idence of a zero-T transition. In addition, Monte-,Carlo 
simulations by Fisher et alB and Reger and YoungQ show 
evidence of a zero-temperature transition, although the 
simulations were not taken down to extremely low tem- 
peratures. 

In this work we perform Monte Carlo simulations of 
the-, two-dimensional gauge glass, using parallel temper- 
ingpEI to go to significantly lower temperatures than was 
possible in earlier work. In particular we are now able 
to cover the temperature range where the claimed spin- 
glass transitionna takes place. We find strong evidence 
that T c = 0. 

In addition, we study the gauge glass in three dimen- 
sions at very low, but finite, temperatures to provide 
a good estimate of the stiffness exponent 9. The mo- 
tivation for this is that earlier estimates, which were 
obtained fp©jp-||g-s©«nd state methods, are inconsistent. 
One groupl!300y finds values consistent with 9 w 0, 



elusions are summarized in Sec. VI 



II. 



MODEL, OBSERVABLES AND 
EQUILIBRATION 



The Hamiltonian of the gauge glass is given by 



H = 



-j 2_, cos ( 

(i,3) 



(1) 



where the sum ranges over nearest neighbors on a square 
lattice in D dimensions of size N = L D and fa represent 
the angles of the XY spins. Periodic boundary conditions 
are applied. The are quenched pandom variables uni- 
formly distributed between [0, 2tt]£3 Because Ay repre- 
sent the line integral of the vector potential between sites 
% and j, we have the constraint that Ay = — A^, In this 
work we set J = 1 . j— . 

Traditionally one uses the Binder ratidl3 to estimate 
the critical temperature T c . In this method one plots 
the ratio of the fourth moment and the second moment 
squared of the spin-glass order parameter as a function 
of temperature for different system sizes L. The crossing 
point of all curves identifies T c . For the gauge glass the 
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Binder ratio cannot exceed unityJH thus the splaying of 
the curves is small and T c is difficult to establish. 

In order to avoid this problem we use a method intro- 
duced by Reger and YoungEj in which one calculates the 
current / defined as the derivative of the free energy F 
with respect to an infinitesimal twist to the boundaries, 



sin 



A-i %-\-x ) 



(2) 



In order to compare with the results of Choi and ParkEl 
we have also calculated the spin-glass susceptibility x SG 
defined by 

Xso = N[(q 2 )U , (4) 

where q is the spin-glass order parameter defined by 



f 



N 

g = -5>xp[i(<^ 



(5) 



In this case the twist is applied along the x direction. 
As the gauge fields Aij are uniformly distributed we ex- 
pect [(7(L))] av = 0, where [• ■ -] av represents an average 
over disorder and (■ ■ •) represents a thermal average. We 
calculate the root-mean-square current 



In 



[{I a (L)){I p (L))] t 



(3) 



Here a, (3 denote two replicas of the system with the same 
disorder introduced to avoid any bias in the calculation 
of the average of the square of the currents. Unlike the 
Binder ratio, the currents are not restricted to converge 
to unity for T 0, thus any potential splaying of the 
data is easily visible. In fact, for T < T c the currents 
scale as / rms ~ L e where 9 is the (positive) stiffness ex- 



ponent, see Refs. p0| , piy22| for a general discussion of the 
stiffness exponent in the context of the "droplet model" 
of spin glasses, and also Ref. ^3| which discusses the stiff- 
ness exponent specifically in the context of the vortex 
glass. 

Note that the stiffness exponent is also often defined 
in terms of the free energy change for a finite twist, by tt 
for example. We are assuming here that any reasonable 
measure of the stiffness will give the same exponent 9, 
whether it be the free energy change due to a finite twist 
or, as here, the derivative with respect to an infinitesimal 
twist. This is consistent with the usual assumptions of 
universality for phenomena on long length scales, that the 
same exponent will be obtained for the relation between, 
say, energy and length scale, for any reasonable definition 
of energy. Presumably a plot of the free energy against 
twist angle O will have quite complicated structure, with 
local regions varying quadratically as a function and 
proportional to L d ~ 2 (where only spin- wave excitations 
are generated) , but with "kinks" where which the system 
changes from one parabola to another due to a vortex 
changing position. Hence, in effect, we are assuming that 
the typical magnitude of the slope of this curve is of the 
same order in L as the range of its variation with 0. 
This requires that the typical spacing between kinks A0 
satisfies L e ~ (A0) 2 i d " 2 , i.e. A0 ~ L-( d - 2 - e )/ 2 . 

For a finite-temperature transition, T c > 0, we expect 
9 > 0, since then the ordered state at T = is "stiff" on 
large scales, and so will presumably resist small thermal 
fluctuations. On the other hand for a zero-temperature 
transition, we expect 9 < 0, because the system will then 
easily break up under the influence of thermal fluctua- 
tions. 



and a, (3 are two replicas of the system with the same 
disorder. This susceptibility di verg es at the transition. 
However, as we shall see in Sec. Ill, its finite-size scaling 
behavior is more complicated than that of I Tms . 

For the simulations, we use the parallel tempering 
Monte Carlo methocoEl as it allows us to study systems 
at lower temperatures than with conventional methods. 
Because the equilibration test fopshort-range spin glasses 
introduced by Katzgraber et al£3 does not work for the 
gauge glass since the disorder is not Gaussian, we test 
equilibration by the traditional technique of requiring 
that different observables are independent of the num- 
ber of Monte Carlo steps N svreep . Figure |l| shows data 
in 2D for the energy E, spin-glass susceptibility x SG , as 
well as the average current squared 7 r 2 ms as a function of 
Monte Carlo steps. One can clearly see that the differ- 
ent observables saturate at the same equilibration time. 
We show data for an intermediate size as it allows us to 
better illustrate the procedure by simulating with much 
longer runs than are necessary to equilibrate. For large 
sizes we equilibrate doubling the number of Monte Carlo 
sweeps between each measurement until the last three 
agree within error bars. 

We also require that the acceptance ratios for the 
global moves which interchange the different tempera- 
tures in the parallel tempering scheme be greater than 
0.3 on average and roughly constant as a function of tem- 
perature. The equilibration tests for the 3D data have 
been described elsewhere.ta 

Tables g and || show the number of samples N samp and 
the number of Monte Carlo sweeps N sweep performed by 
each replica for each lattice size for D = 2 and D = 3, re- 
spectively. In two dimensions the highest temperature is 
1.058, whereas the lowest temperature is 0.13 (for L = 24 
the lowest temperature is 0.20). The number of temper- 
atures Nt is chosen to give satisfactory acceptance ratios 
for the Monte Carlo moves between the temperatures. In 
three dimensions the lowest temperature studied is 0.05 
[to be compared with T c ss 0.45, (Ref. |l|)] whereas the 
highest temperature is 0.947. 

Because the gauge glass has a vector order parame- 
ter symmetry, to speecLup-the simulations we discretize 
the angles of the spinsE3't3 to N,p = 512. To ensure a 
reasonable acceptance ratio for single-spin Monte Carlo 
moves, we pick the proposed new angle for a spin within 
a temperature-dependent acceptance window about the 
current angle. By tuning a numerical prefactor we en- 
sure the acceptance ratios for these local moves are not 
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FIG. 1: Average current squared 7^ ms , energy E and spin- 
glass susceptibility x SG as a function of Monte Carlo sweeps 
A'swccp, that each of the replicas perform, averaged over the 
last half of the sweeps. Note that the different observables 
equilibrate roughly at the same time and appear to be in- 
dependent of the number of sweeps. The data shown is for 
D = 2, L = 6 and T = 0.13, the lowest temperature studied 
in two dimensions. 



TABLE I: Parameters of the simulation in two dimensions. 
A?samp is the number of samples, i.e., sets of disorder realiza- 
tions, iV SWCC p is the total number of sweeps simulated for each 
of the 2Nt replicas for a single sample, and Nt is the number 
of temperatures used in the parallel tempering method. For 
L < 16 the lowest temperature is 0.13, while for L = 24 it is 
0.20. 



L N samp ^Vj sweep 

1 10400 8.0 x 10 4 30~ 

6 10150 8.0 x 10 4 30 

8 8495 2.0 x 10 5 30 

12 6890 8.0 x 10 5 30 

16 2500 2.0 x 10 6 30 

24 2166 2.0 x 10 6 24 



smaller than 0.2 for each system size at the lowest tem- 
perature simulated. 

III. FINITE-SIZE SCALING 

We review briefly some salient aspects of finite-size 
scaling theory, as applied to the present situation, pay- 
ing particular attention to features which are special to a 
T = transition. Consider first the singular part of the 
bulk free energy per spin divided by T, a quantity with 
dimension (length) ~ D , where D represents the space di- 
mension. The basic assumption is that the only relevant 
lengthscale is the correlation length £, and so (3f s ~ (,~ D ■ 



TABLE II: Parameters of the simulation in three dimensions. 
The lowest temperature studied is T = 0.05. 



L 


N samp 


N sweep 


Nt 


3 


10000 


6.0 x 10* 


53 


4 


10000 


2.0 x 10 4 


53 


5 


10000 


6.0 x 10 4 


53 


6 


5000 


2.0 x 10 5 


53 


8 


2000 


1.2 x 10 6 


53 



Consider next a finite-size system where the bulk behav- 
ior is assumed to be modified by a function of L/£, where 
L is the lattice size, i.e., 

ps* = r D f (|) • (6) 

Thus, we can write the scaling of the total singular free 
energy F s (= L D f s ) as 

pF,=p(j) =F{L 1 '»{T-T C )). (7) 

The last expression follows because £ ~ (T — T c )^ lJ and 
we have taken the argument of the function F to the 
power \ jv. 

As J xms is the derivative of F s with respect to the twist 
angle and since the twist angle is dimensionless, the scal- 
ing of i rms is the same as that as of F a , i.e., 

pi Tna =t[L 1 /"(T-T c )]. (8) 

Note that we have been careful to maintain the factor 
of (3 in our analysis. However, if T c > 0, then, since the 
critical region is close to T c , we can replace f3 by 1/T C and 
incorporate this constant factor into the scaling function, 
i.e., 

I rna = l[L 1 ' v {T-T c )\ (T c >0). (9) 

By contrast, if T c — we have to keep the variation in /3 
and so 

7 rms = Ti{L x l»T) 

= L- 1 ^I(L 1 ^T) (T c = 0) , (10) 

where the scaling functions I and / are simply related 
to each other. Equation ( |Io| ) indicates that the T = 
stiffness exponent 9 is negative and equal to —1/v. 
Equation @ shows that if T c > the curves for 7 rms 
for different sizes intersect at the critical point, whereas 
Eq. (|l(]) shows that if T c — 0, then the data decrease 
with increasing size at T = 0. 

According to standard finite-size scaling the spin-glass 
susceptibility, defined in Eq. (Q), behaves as 

X SQ =L^x SG [L 1/v {T-T c )}, (11) 

which means that at criticality it diverges with a power 
law, i.e., 

XsG ~ L 2 -" (T = T c ) . (12) 
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FIG. 2: Root-mean-square current / rms in two dimensions 
as a function of temperature for different system sizes. At all 
temperatures the data decrease with increasing L indicating, 
from Eq. {J9[) , that if T c is finite it must be less than 0.13. 




0.5 1 

L i/v T 

FIG. 3: Scaling plot of the root-mean-square current / m! 
in two dimensions according to the form expected if T c = 
0, Eq. ([lo]). We see acceptable scaling of the data at low 
temperatures. Deviations at higher T are presumably due to 
corrections to scaling. This plot is for 9 = — 1/v = —0.39. 



Equations (|Tl]) and (|l2|) are valid whether the transi- 
tion is at T c — or T c > 0. The power law prefactor 
in Eq. ( |ll"| ) with an unknown exponent complicates the 
analysis of x SG compared with that for / rms which does 
not have such a prefactor if T c > 0, see Eq. (^). Even 
if T c = 0, where a prefactor does arise [see Eq. (juj)], 
the exponent in the prefactor is the same as that in the 
argument of the scaling function, so we do not have an 
extra exponent to determine. 



IV. RESULTS FOR D = 2 

In Fig. ^| we present results for Z rms in two dimensions. 
At all temperatures the data decrease with increasing L 
indicating, from Eq. (^J), that T c must be less than the 
range of temperatures studied. In fact, we shall find it 
impossible to scale the data far I rms for any finite T c , in 
contrast to the results of Kimcl and Choi and Park,u who 
find a finite-temperature transition at T c = 0.22. 

In Fig. U we show a finite-size scaling plot of the data 
in Fig. H according to Eq. (10). Although the data do 
not scale well over the whole range, the data at low-T 
do scale quite well, since the plot shows that the data 
collapse if L X I V T is small for all values of L, and over the 
whole range of L X I V T for the largest sizes. The data in 
Figs. |2| and |^ are therefore consistent with a zero tem- 
perature transition, but with significant corrections to 
scaling at intermediate temperatures. We estimate the 
stiffness exponent to be 



6 = -0.39 ±0.03 



(13) 



The above error bar is estimated by varying 9 slightly 
until the data do not collapse well. This result,is consis- 
tent with recent work of Akino and KosterlitzO who find 
9 = -0.36 ±0.01. 

Figure || shows the current data according to Eq. (||) 
assuming the parameters found by Choi and ParkB: T c = 
0.22, 1/V = 0.88. The scaling is very poor especially 
near the proposed T c . By contrast, if the deviations are 
due to corrections to scaling, they should be smaller near 
T c . Hence our data are incompatible with the claim of 
KimB and Choi and ParkB that T c = 0.22. In fact, we 
are unable to get a reasonable fit to the data for 7 rms 
accor ding to Eq. (|) for any finite T c . 

Next we present our data for the spin-glass suscep- 
tibility x S g ■ Figure |5| is a log-log plot of x SG vs L a ^ 
several low temperatures. According to Eq. (|l^) the 
data should lie on a straight line at T c . However, the 
data in the vicinity of X, = 0.22, the transition tem- 
perature claimed by KimH and Choi and Park,B show a 
strong downward curvature, indicating that this is actu- 
ally above T c . Only around the lowest temperature where 
we have data, T = 0.13, is the curvature small, although 
it still greatly exceeds the error bars. This indicates that 
T c < 0.13, which is compatible with our data for I rms . 

Figure ^ shows a scaling plot according to Eq. (|Tl] ) for 
T c = 0, and 1/v = 0.39, the same parameters found in 
the scaling of J rms , together with r/ = which is expected 
at a zero-temperature transition in two dimensions. The 
data at low temperatures and for the largest sizes scale 
well, but the data away from this range show deviations. 
Allowing 1/v to vary we find 1/v = 0.50 ± 0.03. The 
inset shows data for the optimal value, where only the 
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FIG. 4: Scaling of the current according to Eq. (0) with 
the parameters of Choi and Park (Ref. ^): T c = 0.22, and 
l/u — 0.88. One can clearly see that the scaling is very poor, 
especially in the vicinity of the proposed T c . 
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FIG. 6: Scaling of the spin-glass susceptibility Xsg according 
to Eq. ( |rij ) with T c — 0. The data for large sizes and low 
temperatures collapse with r\ — and l/u — 0.39. The inset 
shows a scaling plot for the optimal value 1/v = 0.50 (and 
?7 = 0). For these values of exponents the collapse extends to 
a larger range of sizes. 
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FIG. 5: Log-log plot of the spin-glass susceptibility Xsg 
as a function of system size L for several different low tem- 
peratures. We clearly see a downward curvature in the data 
indicating that T c < 0.13. 
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FIG. 7: Scaling of the spin-glass susceptibility x S g accord- 



ing to Eq. (|llj) for T c = 0.2S 
suggested by Choi and Park£ 



1/v = 0.88 and r) = 0.30, as 



L = 4 and 6 data are not part of the scaling function 
for all L X I V T. The fact that the best values of 1/v are 
not precisely the same when obtained from x SG and J rms 
presumably indicates that scaling is only valid for fairly 
low temperatures and large sizes, and that, despite our 
working at quite low temperatures, we have only a limited 
range of data which are fully in the scaling regime. 
In Fig. ^ we show the same data as in Fig. ^| scaled 



with the parameters used by Choi and Park in Ref. 
We note that the quality of the data collapse is poor near 
T c , whereas if the deviations were due to corrections to 
scaling, we would expect it to be best in this region. Wc 
therefore do not consider this scaling to be acceptable. 
Furthermore, the scaling of I rms shown in Fig. [| is clearly 
much worse, indicating that I lms is much better able to 
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FIG. 8: Scaling of the spin-glass susceptibility Xsg accord- 
ing to Eq. ( |TT| ) for T c = 0.13, the lowest temperature in our 
simulations, \jv = 0.68 and n = 0.19. The inset shows the 
scaling of the I Tms according to Eq. One can clearly see 
that this data do not scale. 
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FIG. 9: Log-log plot of / rms vs L for different low tem- 
peratures. The data fit well to aL with fitting probabilities 
between 0.57 for the lowest temperature T = 0.050 and 0.35 
for T = 0.176. 



distinguish between a finite T c and T c — than is x SG 
because of its simpler finite-size scaling form. 

This is further illustrated in Fig. |s| which shows a scal- 
ing plot of Xsg m which we took T c to be the lowest 
temperature in our simulations, i.e., T c = 0.13. The best 
fit with this T c has 1/v = 0.68 and r\ = 0.19, where ij 
is determined by the requirement that the data scale as 
well as possible at T c . By eye, the data scale fairly well, 
though the deviations are actually significantly greater 
than the error bars, as we noted in the above discussion 
of the unsealed x S g data at T = 0.13. However, the data 
for Irmsi scale very badly with T c — 0.13 as can be seen 
in the inset of the figure. In fact, as noted above, the 
data for I rms do not scale according to Eq. @ for any 
finite T c . 

The data for x SG enable us to determine the exponent 
r\ as well as \jv. Assuming T c — 0, as obtained from 
the 7 rms data, we find r\ = 0.0 ± 0.1 from a scaling plot 
for x SG . This is consistent with r\ = 0, as expected at a 
T = transition in two dimensions, at least if the ground 
state is not highly degenerate. 



L e when L X / V {T — T c ), the argument of the scaling func- 
tion in Eq. (|^), tends to infinity. The data are shown 
in Fig. ^, in which one observes that the slopes in the 
log-log plot vary somewhat with temperature. Asymp- 
totically, however, one expects the same L 6 behavior for 
all T < T c , so the deviations from a constant slope pre- 
sumably indicate that that the data are not yet at suf- 
ficiently large sizes and low enough temperature to be 
in the asymptotic region. To obtain an estimate of 9 
we therefore perform a linear least-squares fit of In 7 rms 
against In L for each temperature in order to obtain an 
effective stiffness exponent 9 Q g (T) which depends on the 
temperature. Figure |l^ shows that 9 e g(T) can be fitted 
well to a linear form at low temperatures. Extrapolating 
to T — we obtain 

9 = 0.27 ±0.01 (14) 

which is clea rly po sitive and consistent with the results 
ofRefs. |0|5P||. 

VI. CONCLUSIONS 



V. RESULTS FOR D = 3 

The critical region of the three-dimensional gauge glass 
has been investigated in detail by Olson and YoungEl In 
their work they obtain a lower bound for the, stiffn c^ 
exponent of 9 > 0.18. Some previous results— — 
find 9 in the range < 9 < 0.077 whereas othersliJO 
find a much larger value, 0.26 < 9 < 0.31. 

We can estimate 9 from our data for 7 rms since I TZ 



We have shown results from Monte Carlo simulations 
at low temperatures of the two- and three-dimensional 
gauge glass. In two dimensions our data for 7 rms are 
consistent with a zero-temperature transition with 9 = 
—0.39 ±0.03. However, it is necessary to go to quite low 
temperatures and large sizes to see the expected scaling. 
Ourpdata are incompatible with the prediction made by 
KimH and Choi and ParkB that T c = 0.22. The error bar 
we quote is purely statistical and systematic corrections 
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FIG. 10: Effective stiffness exponent 8 a s(T) as a function of 
temperature in three dimensions for low temperatures. The 
data extrapolate linearly to T = 0. 



to scaling can increase this. Indeed, our best estimate 
for \jv = —9 from the Xsg data is 0.50 ± 0.03, which 
differs from the value of 1/v from I tms by more than the 
(statistical) error bars. However, our claim that T c = 
is very robust and is not affected by corrections to finite- 
size scaling. 

In three dimensions we report the first reliable esti- 
mate of the stiffness exponent from finite-temperature 
Monte Carlo simulations. Monte Carlo has the advan- 



tage over ground state methods that it computes the stiff- 
ness directly and is free from the difficulty of determining 
the optimal boundary condition changes, which has been 
quite controversial for the ground state approach. We 
find 9 = 0.27 ± 0.01, which agrees with the results of 
Refs. HlflT5l,|iT^ 
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